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1.   INTRODUCTION 

Ford  and  Fulkerson  (1956)  observed  that  an  easy  way  to  find  the 

value  of  the  maximal  possible  flow  through  an  undirected  source-sink 

* 

planar  network  is  to  construct  the  dual  graph,  assign  the  capacities 

of  the  intersected  primal  arcs  as  lengths  of  the  undirected  arcs  in 
the  dual,  and  then  find  the  shortest  route  through  the  dual.   The 
length  of  this  shortest  route  will  be  equal  to  the  value  of  the  mini- 
mum cut  set  of  the  primal  network  and  hence  is  equal  to  the  value  of 
the  maximal  flow  because  of  the  min-cut  max-flow  theorem. 

Extension  of  this  idea  to  directed  networks  began  with  de 
Ghellinck  (1961) .   He  developed  a  convention  for  constructing  a  dual 
of  a  directed  source-sink  planar  network  but  was  not  concerned  with 
arc  lengths.   In  a  paper  written  in  1969,  Doulliez  and  Rao  (1971) 
provide  the  dual  construction  and  the  arc  length  conventions  for 
directed  capacitated  flow  networks  when  zero  lower  bounds  exist  on 
arc  flows.  This  convention  was  apparently  the  result  of  conversations 
with  de  Ghellinck. 

Because  of  a  need  for  dual  graphs  of  directed  network  in  studies 
of  network  interdiction  models,  the  author  (1970)  independently  devel- 
oped the  same  convention  as  Doulliez  and  Rao  but  went  one  step  further 
and  provided  a  convention  for  handling  problems  where  the  lower  bound 
on  arc  flow  could  take  any  negative  value  or  any  positive  value  not 
exceeding  the  upper  bound  capacity.  This  paper  also  observed  that  for 


Ford  and  Fulkerson  (1956)  defined  source-sink  planar  networks  as  planar 
networks  which  remain  planar  after  an  arc  connecting  the  source  and  sink 
is  added  to  the  network. 


flow  networks  having  positive  lower  bounds  on  arc  flows,  an  infeasible 
primal  problem  is  indicated  by  a  cycle  of  negative  length  in  the  dual 
network.   In  a  later  paper  (1971)  the  author  proved  that  a  cycle  of 
negative  length  was  a  necessary  and  sufficient  condition  for  infeasi- 
bility.   In  addition,  a  proof  was  provided  that  the  value  of  the  minimal 
feasible  flow  could  be  obtained  by  changing  the  sign  on  the  arc  length 
of  the  dual  and  finding  the  longest  route  from  the  destination  back 
to  the  origin. 

During  1970,  M.  F.  Sakarovitch  was  also  working  on  these  same 
problems.   His  research  (1970)  was  motivated  by  some  conversations 
with  de  Ghellinck.   He  also  proved  that  a  cycle  of  negative  length  was 
a  necessary  and  sufficient  condition  for  inf easibility .   In  addition, 
he  proved  that  an  optimal  flow  in  a  primal  arc  of  a  maximal  flow  net- 
work is  equal  to  the  difference  between  the  dual  shortest  distances 
from  the  origin  to  the  two  nodes  at  opposite  ends  of  the  dual  arcs 
intersecting  the  primal  arc.   This  important  property  will  play  a 
major  role  in  the  algorithm  to  be  presented  below  for  non-planar 
directed  networks. 

All  of  the  work  just  summarized  is  restricted  to  those  networks 
which  are  source-sink  planar.   The  main  purpose  of  this  paper  is  to 
extend  these  ideas  to  general  flow  networks. 

In  the  sections  to  follow,  we  will  first  review  the  "state  of 
the  art"  for  source-sink  planar  directed  network  flow  problems.   We 
will  then  turn  to  the  problem  of  non-planar  networks.   A  pseudo  dual 
will  be  defined  and  the  special  form  of  the  dual  shortest  route  problem 
will  be  stated.   The  details  of  the  algorithm  for  solving  this  problem 


will  then  be  presented  followed  by  an  example,  and  the  proofs  which 
validate  the  algorithm.   The  final  section  presents  a  brief  discussion 
of  the  usefulness  of  the  algorithm. 


2.   MAXIMAL  FLOW  IN  A  PLANAR  NETWORK 

The  original  flow  network  will  be  called  the  primal  network. 
A  mesh  of  a  planar  primal  network  is  any  region  surrounded  by  nodes 
and  arcs  but  containing  neither  in  the  plane  on  which  the  network  is 
constructed.   The  region  of  the  plane  completely  surrounding  the 
primal  network  will  be  called  the  external  mesh. 

The  construction  of  the  dual  of  a  source-sink  planar  directed 
network  consists  of  the  following  steps  (McMasters  (1970)): 

1.  Connect  an  artificial  arc  between  the  sink  and  source 
of  the  primal  and  position  it  below  the  network  such  that  it 
crosses  no  other  arcs.  The  resulting  network  will  be  referred 
to  as  the  modified  primal  network. 

2.  Place  a  node  in  each  mesh  of  the  modified  primal  in- 
cluding the  external  mesh.   Let  the  origin  of  the  dual  be  the 
node  in  the  mesh  involving  the  artificial  arc  and  the  destina- 
tion be  the  node  in  the  external  mesh. 

3.  For  each  arc  in  the  primal  (except  the  artificial  arc) 
construct  two  oppositely  directed  arcs  that  intersects  it  and 
join  with  nodes  in  the  meshes  adjacent  to  it. 

A.   Assign  the  value  of  the  upper  bound  capacity  of  the 
primal  arc  as  the  length  of  the  intersecting  dual  arc  having 
the  same  direction  that  the  primal  arc  would  have  if  it  were 
rotated  90  degrees  counterclockwise.   Assign  to  the  oppositely 
directed  dual  arc  a  length  equal  to  the  negative  of  the  lower 
bound  capacity  of  the  primal  arc. 


The  convention  described  in  step  4  is  illustrated  in  Figure  1 

where  the  arc   (i,j)   is  the  primal  arc  having  its  flow  X..   constrained 

ij 

between  two  values  L..   and  M..   such  that  L   £  X   s:  M   .   The 

x3  ij  ij  ij  ij 

adjacent   dual  nodes   are     u     and     v.      Arc      (u,v)      is  assigned  a  length 

equal   to     M  while  arc      (v,u)      is   assigned  a   length  equal   to      -L . . . 

We  assume     M       ^  0     but  place  no   restriction  on     L..      except   that   it 
ij  xj 

cannot   exceed     M.  .  . 
3-1 


fuv=Mjj 


Afu  =  ~Ljj 


Figure  1. 

Figure  2a  is  an  example  of  a  flow  network  having  node  1  as  its  source 

and  node  4  as  its  sink.   The  numbers  on  each  are  represented  L... 

ij 

M. . .   Figure  2b  shows  the  initial  phase  of  the  construction  of  the 
dual.   The  dual  nodes  are  designated  as  A,   B,   C  and  D.   Each 
dashed  arc  corresponds  to  the  location  of  each  pair  of  dual  arcs.   Fig- 
ure 2c  is  the  completed  dual  with  its  associated  arc  lengths. 

Ford  and  Fulkerson  (1962)  defined  the  generalized  value  of  a 
primal  cut  set  when  lower  bounds  are  positive  by  equation  (1)  where 


Figure  2a. 


Figure  2b 


Figure  2c . 

the  set   S   contains  the  primal  source  node  and  the  set   S   contains 
the  primal  sink  node. 


V(S,S)  =  I     M..  -  I      L...  (1) 

ieS   ^    ieS   1J 
jeS       jeS 

They  then  extended  the  min-cut  max-flow  theorem  to  include 

this  definition  of  the  value  of  a  cut  set;  namely,  if  a  feasible  flow 

can  be  found  for  a  network  with  nonnegative  lower  bounds  on  arc  flows 


then  the  maximal  possible  flow  through  the  network  is  equal  to  the 
minimal  generalized  value  of  all  possible  cut  sets  of  the   (S,S)   type. 
Later  Fulkerson  (1962)  dropped  the  restriction  that  the  lower  bound 
be  nonnegative. 

The  relationship  between  the  convention  of  Figure  1  and  the 
generalized  value  of  a  primal  cut  set  is  described  by  the  following 
theorem: 

Theorem  1 :  Any  dual,  chain  without  cycles  directed  fatiom 
the,  dual  houKce  to  the,  dual  6ink  hcu  a  length  equal  to 
the  QontnaJLLznd  value.  o£  the  primal  cut  6et  intersected 
by  the  axch  o{  the  chain. 

Proof:   Each  path  without  circuits  through  the  dual  of  a  source- 
sink  planar  undirected  graph  intersects  a  cut  set  of  the  primal  network 
such  that  the  primal  source  and  sink  are  separated  by  the  cut  and 
there  is  a  separate  path  corresponding  to  each  cut  (Ford  and  Fulkerson 
(1956),  Whitney  (1933)).   This  is  still  true  if  the  primal  and  dual 
arcs  are  replaced  by  directed  arcs.   Further,  because  the  dual  graph 
under  the  Figure  1  convention  has  two  oppositely  directed  arcs  between 
every  pair  of  nodes  in  adjacent  primal  meshes  it  is  also  true  that  a 
separate  chain  exists  directed  from  the  dual  source  to  the  dual  sink 
corresponding  to  each  primal  cut. 

Suppose  we  select  an  arbitrary  cut  set  of  the  primal  network 
such  that  the  source  is  in  a  set   S   and  the  sink  is  in  set   S.   Any 
such  cut  set  can  be  illustrated  as  shown  in  Figure  3  for  planar  graphs. 


Notice  that  this  set  contains  both  arcs  directed  from  S   to  S  and 
from  S   to  S.   Let  the  chain  directed  from  the  dual  source  to  the 
dual  sink  which  intersects  the  cut  set  arcs  be  a-b-c-d-e-f 


DUAL  SINK 


Figure  3. 


Using  the  convention  of  Figure  1  we  assign  a  length  of  -L 
to  dual  arc   (a,b)  which  crosses  primal  arc  t,   a  length  of  M   to 
arc   (b,c)   which  crosses  arc   s,   a  length  of  M   to  arc   (c,d), 
a  length  of  -L   to  arc   (d,e),  and  a  length  of  M   to  arc   (e,f). 
The  chain  directed  from  a  to   f  would  therefore  have  a  length  of 


M  +  M  +M  -L   -L 
s    r    p    t    q 


which  we  see  immediately  corresponds  to  the  form  of  the  right-hand 
side  of  (1)  for  this  arbitrary  cut  set  (the  arcs  p,   r  and  s  have 
ieS  and  jeS  while  arcs   q  and   t  have  ieS  and  jeS).   Because 
our  cut  was  arbitrary  we  realize  that  if  we  analyze  any  cut  set  of  the 
(S,S)   type  in  this  manner  we  will  get  the  length  for  its  dual  chain 


being  described  by  this  general  form.   Therefore,  each  chain  without 
cycles  from  the  dual  source  to  the  dual  sink  will  have  a  length  equal 
to  the  generalized  value  of  the  primal  cut  set  it  intersects. 

As  a  consequence  of  Theorem  1  and  the  generalization  of  the 
min-cut  max-flow  theorem  we  get  the  following  corollary: 

Coh.oHaAy  7:     The.  value.  o&  maximal  falovo  In  a  ^exulble. 
dunuLtud  Aou/L(ie.-6Znk  planan.  network  u>  equal  to  the. 
length,  o&  the.  i>hon£eAt  noute.  tfoiouah  the.  dual  network. 

Recalling  Figure  2a  we  know  that  the  value  of  the  primal  cut 
set  consisting  of  arcs   (1,2),   (2,3),   and   (3,4)   is 

M12  +  M34  "  L23  =  5  +  3  "  X   =  7*  (2) 

The  route  through  the  dual  network  which  corresponds  to  this  cut  set 

is  the  chain   [ (A,C) , (C,B) , (B,D) ]   which  has  a  length  given  exactly 

by  equation  (2) .   This  chain  happens  to  be  a  minimum  route  of  Figure 

2c  so  that  the  maximum  possible  flow  through  the  primal  network  is, 

in  fact,  7  units.  As  we  will  show  below,  the  fact  that  we  found  a 

shortest  route  of  finite  length  means  that  the  primal  problem  is  feasible, 

The  extension  of  the  min-cut  max-flow  theorem  involving  the  use 
of  the  generalized  value  of  a  cut  set  as  given  by  (1)  is  contingent 
on  the  existence  of  a  feasible  flow.   Feasibility  requires  that  a  set 
of  arc  flows  can  be  found  which  satisfy  the  flow  conservation  equations 
at  each  node  as  well  as  the  upper  and  lower  bounds  on  individual  arc 
flows.   When  the  lower  bound  is  nonpositive  then  we  know  that  a  flow 
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problem  will  always  be  feasible  since  X. .  »  0   is  a  feasible  solution 
under  the  M. .  ^  0  assumption.   If,  however,  the  lower  bounds  on  some 
arcs  are  positive  then  we  cannot  claim  such  problems  will  be  always 
feasible.   An  example  of  an  infeasible  problem  can  be  obtained  by 
changing  the  lower  bound  on  arc   (1,3)   in  Figure  2a  to  a  value  of 
3. 

The  answer  to  the  feasibility  question  is  provided  by  Theorem  2, 


Theorem   2:  A  nccc66aAy  and  Au^tctent  condition  {on.  an 
tn{cai>lblt  {low  problem  abbociated  with  a  directed 
6ouA.cc-6infz  ptanaA  network  -U  tkout  a  cycle,  ofi  negative 
length  cxAAtA  in  the  dual  netitionk. 

The  proof  of  Theorem  2  depends  upon  the  circulation  theorem  of 
Hoffman  (1960)  which  is  addressed  to  networks  having  a  sourceless- 
sinkless  structure.   Such  a  network  is  easily  created  from  any  flow 
network  by  connecting  a  directed  arc  from  the  sink  to  the  source  and 
specifying  the  appropriate  upper  and  lower  bounds  on  its  arc  flow. 

The  circulation  theorem  states  that  a  network  flow  problem  will 
be  feasible  if  and  only  if 


I     \i    *  L  Lij  (3) 

ieY   J    ieY   J 
jeY       j£Y 

for  all  cut  sets   (Y,Y)   of  the  network  after  it  has  been  converted 

to  the  sourceless-sinkless  form. 


Suppose  we  define  the  value  of  the  cut  set   (Y,Y)   as 

V(Y,Y)  -  I     M   -  £  L   .  (4) 

ieY   J    ieY   J 
jeY       jeY 

If  we  do  this  then  the  feasibility  condition  given  by  (3)  is  equiva- 
lent to  V(Y,Y)  ^  0  and  inf easibility  corresponds  to  some  V(Y,Y)  <  0, 

An  interesting  correspondence  between  dual  cycles  and  V(Y,Y) 
is  described  by  Theorem  3.   We  will  prove  this  theorem  before  proving 
Theorem  2  since  it  will  be  useful  in  the  proof  of  the  latter. 

Th.zon.2jm   3;  TkeAe.  -c6  a  i>zpcvioutz  dual  cycle  tnte/u, activiQ 
oaqj>  o{  each  pHAjnal  cut  62X     (Y,Y).  ThJj>  cycle.  hcu>  a 
length  zquml  to     V(Y,Y). 

Proof:   Any  circuit  in  the  dual  corresponds  to  a  cut  set   (Y,Y) 
of  the  primal  (Whitney  (1933)).   Because  of  our  convention  of  having 
two  oppositely  directed  dual  arcs  intersect  each  primal  arc  we  know 
that  two  cycles  exist  (oppositely  directed)  in  the  dual  corresponding 
to  each  cut  set  of  the   (Y,Y)   type  in  the  primal.   Suppose  the  set 
Y   is  circumscribed  by  these  dual  cycles.   The  convention  of  Figure  1 
will  then  result  in  the  cycle  directed  clockwise  having  a  length 
corresponding  to  V(Y,Y).   If  Y   is  the  set  circumscribed  by  the 
cycles  then  the  counterclockwise  cycle  has  a  length  V(Y,Y). 

Figure  4  is  an  example  of  an  arbitrary  cut  set  of  the  (Y,Y) 
type  with  its  associated  dual  cycles.  Primal  arcs  e,   f,  h  and  j 
are  directed  from  Y   to  Y;   arcs  g,   i,   k  and  I     are  directed 
from  Y   to  Y.   The  convention  of  Figure  1  will  result  in  the  cycle 
p-q-r-s-t-u-v-w-p  having  a  length  of 


12 


M     +M.+M+M     -L     -  L-  -  X,     -L 
g  i         Tc  e  f  h  j 


which  is   the  right  hand   side  of    (4) . 


Figure  4 . 

We  will  now  prove  Theorem  2.  Suppose  that  the  primal  problem 
is  infeasible.  From  the  circulation  theorem  we  know  that  there  will 
be  some  cut  set   (Y,Y)   of  the  sourceless-sinkless  primal  which  has 


I  M. .  <  I     L 
ieY   1J    ieY   1J 


(5) 


JeY 


JeY 


We  can  rewrite  (5)  in  the  following  form: 


I  M..  -  l_   L   <  0.  (6) 

ieY   1J    ieY  13 
jeY       jeY 

From  (4)  we  know  that  the  left-hand  side  of  (6)  is  the  value  of  the 

cut  set   (Y,Y).   It  is  also  the  length  of  a  cycle  in  the  dual  because 

of  Theorem  3.   Therefore,  an  infeasible  primal  problem  results  in  a 

cycle  of  negative  length  in  the  dual. 


13 


To  show  sufficiency,  we  consider  a  cycle  of  negative  length  in 
the  dual.   Let  its  length  be  represented  by  D.   Because  of  Theorem  3, 

D"I  M   "  L  L   =  V(Y,Y), 
ieY   J    ieY   J 
jeY       jeY 

and  D  being  negative  means  that  the  value  of  some  cut  set   (Y,Y)   in 

the  primal  is  negative.   For  a  cut  set   (Y,Y)   to  have  a  negative  value 

(5)  must  be  satisfied  and  the  primal  will  be  infeasible. 
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3.      MINIMAL   FLOW   IN  A  PLANAR  NETWORK 

It   is   sometimes  useful   to  also  know  the  minimal   feasible  flow 
for  a  network  having  positive  lower  bounds  on  arc   flows.      The  following 
theorem  provides  an  easy  way  of  determining   such  a  flow. 

Theorem  4:     l&  the  oac  length*  ofi  the  dual  o&  a  feasible 
pnMnal  filou)  network  have  thoJji  signh  sieveteed  then  the 
longest  chain  ctiAected  ^lorn  the  dual  destination  back 
to  the  dual  AouAce  loHZ  have  a  length  equal  to  the  value 
o&  the  minimum  feasible  filau)  through  the  netwoh-k. 

The  proof  of   this   theorem  is   based  on  the  min-flow  equivalent 
of   the  max-flow  min-cut   theorem   (Ford  and  Fulkerson   (1962));    that   is, 


min  Q  =  max       \     L        -     \     M        ,  (7) 

ieS        J        ieS        J 
jeS  jeS 

where     Q     represents   the   total   flow  leaving   the  network. 

The  combination  of   our  convention  of   allocating  dual   arc  lengths 

and   the  sign  reversal  just   described  will   cause  each   chain  from  the 

destination  back  to   the  origin  in  the  dual  network  to  have  a  length 

equal   to 

I      L..  -  I     M... 
ieS   1J    ieS   1J 
jeS       jeS 

The  proof  of  this  property  parallels  that  of  Theorem  1.   Therefore, 

we  get  the  right-hand  side  of  (7)  when  we  find  the  longest  chain  back 

through  the  dual. 
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4.   OPTIMAL  ARC  FLOWS  IN  PLANAR  NETWORKS 

The  final  aspect  of  the  source-sink  planar  network,  analysis  is 
the  determination  of  the  optimal  arcs  flows  of  the  primal  network. 
These  flows  are  available  from  the  solution  to  the  dual  network  shortest 
route  or  longest  route  problem  depending  on  whether  we  want  maximal  or 
minimal  flow  through  the  network. 

Consider  the  maximal  flow  problem.   Let  A  be  the  set  of  all 
dual  arcs  associated  with  a  given  primal  flow  network.   Then  the  dual 
network  shortest  route  problem  can  be  stated  as  follows: 


minimize    T  Y,  I..    U.., 
L   Lk     ij   ij 

(-1  (i  =  D, 

N 
subject   to        I      [U..-U..]=<      0  (i  =  2,3, . . . ,N-1) ,  (8) 

3=1        J1  1J  I 

\+l  (i  =  N), 


and     U. .    ^  0  ((i,j)eA). 

If  we  think  of  U..   and  I,        as  representing  respectively  a 
flow  and  a  unit  flow  cost  through  a  dual  arc   (i,j)   then  we  realize 
that  (8)  corresponds  to  a  minimal  cost  network  flow  problem  for  a 
total  flow  out  of  one  unit.   The  optimal  solution  will  yield  a  chain 
flow  of  one  unit  over  the  cheapest  route  from  the  origin  to  the  desti- 
nation of  the  dual  network.   This  route  obviously  corresponds  to  the 
shortest  route  through  the  dual  network.   The  dual  linear  programming 
problem  associated  with  (8)  can  be  written  as 
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maximize    VM  -  Vn , 
N    1 


subject  to  V.  -  V.  J£  £..      ((i,j)eA). 


(9) 


If  we  apply  any  of  the  current  algorithms  (see  Dreyfus  (1969)) 

for  determining  the  shortest  routes  from  node  1  to  every  other  node 

we  will  obtain  a  number  f   associated  with  the  length  of  the  shortest 

route  from  node  1   to  node  i.   If  we  let  V.  =  f .   then  we  will  have 

l    l 

an  optimal  solution  to  (9) . 

The  relationship  between  the  V.   values  and  the  optimal  arc 
flows  in  the  original  primal  network  is  described  by  Theorem  5  (Sakaro- 
vitch  (1970)). 

Tkz.oK.om   5:  ConiideA  a  primal  one,     (a, 3)  diAecXed  ^Kom 
node,    a  to  nods.     3  uikich  kcu,  a  filou)    x  .  Let  the 

01  p 

inteueeting  dual  ancA  be     (i,j)     and     (j,i)     with  the 

length  o{,     (i,j)     being    M        and  the  length  o&     (j,i) 

ap 

being     -L     .     I A     X        i>atLb{ieA> 

otp  Otp 

X   D   =  V.    -  V.  (10) 

a3  j  i 

faon.  eve/iy  one     (a, 3)     then  an  optimal  {eahible  solution 
haii  been  faound  to  the  pnJjnal  maximal  &low  problem. 


Proof:      Because     V.    -  V.      is   a   feasible  solution  to    (9)   we 

j  i 

know  that 

L   „  S  V.    -  V,   j;  M  „ 
ap1  J  i  ag 

under  our  convention  for  assigning  dual  arc   lengths.      Thus 
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ag    ag    ag 

and  therefore  X    is  a  feasible  flow  through  arc   (a, 3). 

To  check  for  flow  conservation  at  the  primal  nodes  we  first 
introduce  a  new  primal  arc  (n,l)  connecting  the  sink  to  the  source 
and  set  its  upper  and  lower  bounds  as  well  as  the  arc  flow  from  the 

sink  to  the  source  at  the  value  of  V„  -  V, .   We  then  have  a  source- 

N    1 

less-sinkless  structure  and  the  flow  conservation  equations  can  be 
written  as 


I     l"X  a    ~   Xa  1  =  °  (11) 


3 

for  a  =  1,2, .. . ,n. 

If  we  next  add  dual  arcs  to  intersect  the  primal   (n,l) 

arc  we  will  create  a  dual  network  which  has  a  cycle  circumscribing 

every  primal  node.   If  we  assign  arc  lengths  to  these  arcs  according 

to  our  convention,  we  will  be  adding  two  more  constraints  to  (9) . 

However,  the  previous  optimal  solution  to  (9)  is  also  optimal  for 

this  modification. 

If  we  sum  the  V.  -  V.   differences  over  all  arcs  of  any  of 
J    i 

these  dual  cycles  we  will  get  zero  since  each  V.   appears  twice  in 
the  summation,  once  with  a  plus  sign  and  once  with  a  minus  sign. 
Therefore,   X  „  =  V.  -  V.   results  in  the  sum  of  the  flows  entering 
and  leaving  a  primal  node  being  equal  to  zero  and  (11)  is  satisfied. 

Finally,  the  optimal  solution  to  (9)  will  result  in  V.  -  V. 
=  Z..      for  all  arcs  on  the  shortest  route  between  nodes   1  and  N 
of  the  dual  and  equation  (10)  will  result  in  the  flows  in  the  arcs 
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of  the  primal  minimum  cut  set   (S,S)   being  M    if  aeS   and   geS 

ap 

and  L    if  aeS   and   $eS .   No  flow  augmenting  paths  will  therefore 
ap 

exist  between  the  primal  source  and  sink  and  hence  the  existing  flows 
are  optimal.   This  completes  the  proof  of  Theorem  5. 

The  reader  will  recall  from  Theorem  4  that  we  also  can  get 
the  minimal  flow  if  a  problem  is  feasible  by  changing  the  signs  on  the 
arc  lengths  of  the  dual  and  finding  the  longest  route  from  the  sink 
back  to  the  source.   The  associated  longest  route  linear  programming 
problem  can  be  stated  as 


maximize    T  Y ,  I .  .   U .  . , 
L  LA     ij  ij 

N 
subject  to   I      [U   -U ..]={   0     (i  =  2,3, . . . ,N-1) ,     (12) 
j=l    J    J 


and  U . .  £  0 . 

The  corresponding  dual  linear  programming  problem  is 

minimize    W„  -  W. , 
N    1 

subject  to  W.  -  W  >  I..  ((i,j)eA).  (13) 

If  we  determine  the  longest  route  from  each  node  to  the 
origin,  we  will  have  an  optimal  feasible  solution  to  (13) .   The  value 
of  W  -  W   is  minimal  feasible  flow  through  the  network  (a  negative 
value  means  flow  from  the  sink  back  to  the  source) .  The  optimal 
feasible  solution  to  the  original  primal  minimal  flow  problem  can  then 
be  obtained  from  the  results  of  Theorem  6.   The  proof  of  Theorem  6 
parallels  that  of  Theorem  5. 


+1 

(i  =  N), 

0 

(i  =   2,3,...,N-1), 

-1 

(i  =   1), 
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TkcotLem  6:     ZonbidoJi  a  pnJjnoJL  axe     (a, 3)     dJUMLcX<L&  ^Kom 

nodz    a     to  nodo,     g     u)kick  kcu>  a  Alow    X  „.     LoX  the. 

u  ag 

tnteAA toting  dual  a/icA  be.     (i,j)     and     (j,i)     tucth  the. 

lungtk  oh     (i,j)     6ei.ng     -M  0     and  -t/ie  length  ol     (j,i) 

ap 

beXng     L     .     Itf     x        4<LtcA^ce6 

x  0  =  w.  -  w. 
a&         j         i 

faon  zveAy  an.c     (a, 3)     #ien  an  optimal.  heoAtble.  solution 
haA  bzm  hound  to  thz  pulmal  mlntmaji  ^lov)  problem. 
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5.   AN  EXAMPLE 

The  results  of  Theorems  1  through  6  will  be  illustrated  by 
the  following  example.   Figure  5a  is  a  network  having  node  1  as  its 
source  and  node  2  as  its  sink.   The  numbers  on  the  arcs  represent  the 
lower  and  upper  bounds  on  arc  flow.   The  dual  construction  is  shown 
by  the  nodes  A,   B   and  C  and  the  dashed  arcs  connecting  them. 
The  resulting  dual  network  for  the  maximal  flow  problem  is  shown  in 
Figure  5b.   The  maximal  flow  through  the  primal  is  4   corresponding 
to  V  -  V  .   The  dual  shortest  route  is  the  chain  of  arcs  in  Figure 
5b  having  the  hash  marks.   The  primal  minimum  cut  set  is  therefore 
{(1,2), (2, 3)}. 


Vo 


Vc-4 


(b) 


Figure  5. 


The  flow  through  arc   (1,2)   of  the  primal  is   6  as  given  by 

the  difference  V  -  V   since  the  dual  arcs   (B,C)   intersect  it. 
C     o 

The  flow  through  arcs   (2,3)   and   (3,1)   are  both  2   corresponding 

to  the  difference  V.  -  V_ . 

A    B 
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If   L„1   is  changed  to  a  value  of   4   then  the  top   (A,B) 
arc  of  Figure  5b  has  a  length  of   -4.   A  cycle  of  negative  length 
will  exist  consisting  of  the  top   (A,B)   arc  and  the  bottom   (B,A) 
arc.   The  problem  is  infeasible  according  to  Theorem  2.   Inspection 
of  Figure  5a  under  this  change  shows  that  an  inf easibility  condition 
has  been  created  by  arc   (3,1)   having  a  lower  bound  exceeding  the 
capacity  of  arc   (2,3). 

Figure  6  is  the  dual  network  for  the  minimal  flow  problem. 


wc=i 


Figure  6. 


The  minimum  feasible  flow  from  the  source  to  the  sink  of  the  primal 
is  W  -  W  =  1.   The  arc  flows  creating  this  are  X  „  =  W£  -  Wfi  =  4 


aad  X23  -  X   =  WA  -  WB  =  3. 
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6.   MAXIMAL  FLOW  IN  NON-PLANAR  NETWORKS 


We  would  like  to  make  use  of  the  ideas  developed  in  the  preceding 
section  for  all  networks  and  not  just  those  which  are  source-sink  planar. 
The  purpose  of  this  section  is  to  present  a  way  of  approaching  such 
problems. 

Suppose  we  construct  a  two-dimensional  representation  of  a  non- 
planar  network  such  that  all  arcs  are  straight  lines.   If  we  do  this, 
then  no  two  arcs  will  intersect  more  than  once  at  points  which  are  not 
nodes.   If  we  designated  these  extra  arc  intersection  points  as  pseudo 
nodes  we  will  have  created  a  network  which  is  source-sink  planar  and 
we  can  construct  its  dual  network  using  the  procedure  and  convention 
of  section  2. 

Unfortunately,  we  have  no  guarantee  that  the  shortest  route 
through  this  dual  network  will  give  the  maximal  feasible  flow  and  that 
the  associated  optimal  feasible  arc  flows  can  be  determined  from 
Theorem  5.   Consider  Figure  7.   Suppose  arc   (a, 3)   and   (y><5)   intersect 
at  pseudo  node  <J>. 


0 


©■ 


0— 


©• 


PSEUDO  NODE 

-J-O 


<s> 


■® 


© 


Figure   7. 
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If  we  hope  to  use  the  results  of  Theorem  5  to  obtain  the 

optimal  feasible  flows  in  these  arcs,  we  must  require  that  the 

values  V  ,   V.,   V  ,   and  V   in  (9)  for  the  associated  dual  nodes 
i    j     k.         m 


satisfy 


V.  -  V.  =  V.  -  V  , 
i    k    j    m' 


V.  -  V.  =  V  -  V.  . 
J     i    m    k 


(14) 


Condition  (14)  is  needed  because  an  optimal  feasible  solution  to  the 
original  primal  maximal  flow  problem  must  have 

XY<}>  =  V' 

since  arcs   (a,g)   and   (y><5)   do  not  intersect  in  the  original  primal 

network.   Therefore,  for  Theorem  5  to  be  useful,  condition  (14)  must 

be  added  to  the  constraints  of  (9) . 

We  note  in  passing  that  (14)  can  be  consolidated  into  only  one 

equation  which  is 

V.  +  V,  -  VJ  -  V  =0.  (15) 

J    k    i    m 

Let  the  node  group  surrounding  a  primal  pseudo  node  <j>  be 
represented  by   (i,j,k,m)  .   Let   $  be  the  set  of  all  pseudo  nodes 
in  a  given  primal  network.   The  incorporation  of  (15)  into  (9)  then 
results  in  the  following  problem: 
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minimize    V.T  -  V.. 
N    1 

subject  to  V  -  V±  £  I  ((i,j)eA),         (16) 

Vi  +  \  "  Vi    "  Vm  =  °     ((i,j,k,m)  ;   *le«) 
j     k.     i    m  <p 

Obviously,  the  optimal  solution  to  (16)  cannot  always  be  deter- 
mined using  an  algorithm  for  finding  the  shortest  route  from  node  1 
to  each  node  i.   The  following  algorithm  is,  therefore,  suggested  as 
a  way  of  solving  (16) . 

Algorithm: 

1.  Compute  the  shortest  route  lengths  from  the  dual  origin  to  all 
nodes.   If,  in  the  process,  a  cycle  of  negative  length  is  detected, 
terminate  the  algorithm.   The  primal  flow  problem  is  infeasible. 

If  no  cycle  of  negative  length  exists  then  compute  the  value  of 
V.  +  V,  -  V.  -  V   for   (i,i,k,m),  associated  with  each  pseudo 
node  in  $   and  go  to  step  2. 

2.  Examine  the  set  of  V.  +  V.  -  V.  -  V   values  obtained  in  step  1. 

j    k    l    m 

(a)  If  V.+V,  -V.  -V  =0   for  all  pseudo  nodes  in  $   then 

j    k    l    m 

terminate.   The  optimal  solution  to  the  primal  network 
problem  can  be  obtained  by  applying  Theorem  5  to  the  planar 
representation  of  that  network. 

(b)  If  V.  +  V  -  V.  -  V  ?  0     for  any  set   (i,j,k,m)    then 
compute  a  set  of  numbers  Y.   (i  =  1,2,...,N)   where 


YN  '  VN  Md 


Y .  =  max  (Y . -%  . . ) 
1         J   XJ 


for   i  =  1,2,...,N-1.   Go  to  step  3. 
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3.   For  each  set   (i,j,k,m)a   for  which  V.+V,-V,-V  *  0 

9  j    k    i    m 

compute  the  following  numbers : 

*2(*)  =max{-VYWV- 

If  £,(9)  <  £«(<}))   for  one  or  more  sets  then  go  to  step  4; 
otherwise  go  to  step  5. 


4.   Determine  that  set   (i,j,k,m)~   corresponding  to 

9 

£2(<J>)  -  ^(9)  =  maxU2(9)  "  ^(9)  I  ^(9)  <  A2^^ 


and  change  a .  .  ,   £ .  . ,   £,    and  £  .   of  the  arcs  asso- 
xj    j  x    km        mk 

ciated  with  that  set  to 


'ji  ■  •*  ■  ^w- 

Return  to  step  1. 

5.   Select  any  set   (i,j,k,m)   for  which  £,(9)   and  i^) 

values  have  been  computed  and  change  I . . ,  I . . ,      I,         and 

Xj    j  x    km 

I   .   of  the  arcs  associated  with  that  set  to 
mk 

l\.   =  A»  =   I, (9) 
xj    km    1  T 

jx    mk     1 
Return  to  step  1. 
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7.   AN  EXAMPLE 
To  illustrate  the  algorithm  we  will  solve  for  max  Q   in  the 
example  shown  in  Figure  8.   This  network,  while  being  planar,  is  not 
source-sink  planar.   Figure  8  is,  in  fact,  the  two  dimensional  repre- 
sentation.  We  see  that  all  arcs  are  straight  lines  and  one  pseudo 

node  has  been  created.   The  numbers  on  the  arcs  are  the  M. .   values: 

ij 

L..  =  0   for  all  arcs. 

The  dual  network  is  shown  in  Figure  9.   The  V.   and  Y.   values 

are  shown  next  to  their  respective  nodes.   We  see  immediately  that 

V..  +  V  -  V,  -  V  /  0  so  the  Y.   values  are  needed.   Next  we  obtain 
f    e    d    g  i 

2,  (<(>)  =  5  and   £?(<}>)  =  8.   Therefore,  we  change  2,,f   and  I  to 

5,   2,,,   and  I  to   -5  and  recompute  the  shortest  route, 

fd        ge 

The  new  set  of  V.   values  are  shown  in  Figure  10.  We  now  have 

V,.  +  V  -  V,  -  V  =0   so  we  are  done.   The  maximal  possible  flow 
f    e    d    g 

through  the  primal  network  is   7  as  indicated  by  the  difference 

V,  -  V  .   The  optimal  solution  to  the  maximal  flow  problem  associated 
h    a 

with  the  primal  network  is  provided  by  Figure  11.  The  numbers  on  the 
arcs  represent  X../M...   The  minimum  cut  set  is   { (4 ,5) ,  (4,2)  ,  (1,2)  , 
(3, 2), (5, 3)}. 
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SOURCE 


,®- 


SINK 


N>®~ 


Figure  8, 


V10 


Figure   9. 
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Vc«IO 


vfl-o 


Vu«7 


Figure  10. 


SOURCE 


SINK 


Figure   11. 
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8.   PROOF  OF  THE  ALGORITHM 


When  V.  +  V.  -  VJ  -  V  ^0   for  some  set   (i.j.k.mV   we 
jkim  '"**    '      y 

must  seek  new  values,   Vj,   Vj,   V,'  ,   V,  which  will  satisfy  the 

1    j    k    m 

constraints  of  (16) .   In  an  attempt  to  avoid  reducing  the  value  of 

V  -  V   and,  therefore,  the  value  of  max  Q  we  first  compute  the 

Y.   values.   The  differences  YXT  -  Y .   (i  =  1,2,...,N-1)   are  the 

x  N    l 

lengths  of  the  shortest  routes  from  each  of  the  nodes  to  node  N. 

Now 

VN  "  Vi  *  YN  "  V  (17) 

since  V  -  V.  will  be  less  than  or  equal  to  the  length  of  any 

chain  directed  from  node  i  to  node  N   (equality  holds  only  if 

node  i   is  on  the  shortest  route  from  node  1   to  node  N) .   Because 

Y  =  V  ,   (17)  reduces  to 

Y.  <;  V.  . 
l    l 

Therefore,  if  we  select  V!   satisfying  the  constraints  of  (16)  such 


that 


Y.  £  V!  £  V,  (18) 

l    l    i 


then  we  are  assured  that  V  -  V^^  will  not  change. 

Combining  (18)  with  the  inequality  constraints  of  (16)  results 
in  the  following  ranges  for  V'  -  V|  and  V^  -  V^: 

max{-*..,Y.-V.}  £V|  -  V!  £  min{ I      ,V  -Y  } , 


(19) 

<  n,in(P.   .V  -Y.  }. 

mk'  m  k     m    k 


■"<-*.* -w  * v;  - v;  *  min(*km'W- 
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Let  us  assume  the  ranges  in  (19)  have  some  values  in  common. 

We  next  combine  the  inequalities  of  (19)  with  the  equality  constraints 

of  (16) .   The  result  is  (20) .   Notice  that  I   .       and  I,         have  been 

mk        km 

omitted.   This  is  due  to  the  fact  that  I..    -  %    .   and  I . .  =  %   . 

xj    ml  ji    mk 

for  every   (i, j ,k,m)    as  a  consequence  of  the  dual  construction. 


max{-£.  .  ,Y.-V.  ,Y  -V.  }  £  V!  -  V!  =  V*  -  V.'  £  minU  .  .,V.-Y.,V  -Y.  }.    (20) 
ji'ji'mk     j     i    m    k        lj'ji'mk 


We  realize  that  (20)  reduces  to  the  following  form: 


i2c*>  svj  -v^r-v^^), 

when  we  use  the  definitions  of  £..(<{>)   and  l0(<\>)      from  step  2  of  the 
algorithm. 

Step  4  of  the  algorithm  handles  the  problems  where  the  ranges 
of  (19)  have  values  in  common  for  all  $   which  violate  (15) .   It 

changes  the  arc  lengths  of  arcs   (i,j)   and   (k,m)   to  I    (<f>)   and  the 

arc  lengths  of   (j,i)   and   (m,k)   to  -I-  ($).   This  forces  the  condi- 
tions 

V   -  V.  =  V.  -  V.  =  I.  (<j,), 
m    k    j     i    i^t" 

V.  -  V,  =  V.  -  V  , 
x    k    j    m' 

while  satisfying  (18).   Note  that  one  node  of  each  pair,   (i,j)   and 

(k,m) ,   will  have  V'  =  V.   after  the  new  computations.   Finally, 

V'  -  V'  =  VXT  -  V,   so  that  the  value  of  maximal  flow  has  not  been 
N    1    N    1 

affected  by  the  new  computations. 
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If  the  ranges  of  (19)  have  no  values  in  common  for  some  set 

(i,j,k,m)    then  I    ($)  <  fc.(<j>)   and  the  current  V  -  V,   value 
<p  1       I  N    1 

represents  an  infeasible  value  for  max  Q.   Because  V   must  be 

N 

reduced  by  the  amount 

£2(«(.)  -  J^CJ)  =  tnax{&2(<J>)  -  l±W    |  j^Oj,)  <  l2($)} 

before  all  sets   (i,j,k,Jl)    can  meet  the  constraints  of  (16),  step  3 
adjusts  the  appropriate  arc  lengths  of  the  set  or  sets  corresponding 
to  <$>      to   £,(<J>)   and  -I    (<}>)   to  provide  this  reduction  in  V 
immediately.   No  further  use  of  step  3  is  necessary  after  this  reduc- 
tion.  It  is  possible,  however,  that  adjustments  of  the  step  4  variety 
will  still  be  needed  before  an  optimal  solution  is  reached. 

Changing  the  arc  lengths  of  the  arcs   (i,j),   (k,m)   and   (j,i), 
(m,k)   to  i.  (<J>)   and   -&,  (<j>)   values  or  to   £,(9)   and  -£,(<)>)   values, 
respectively,  results  in  the  optimal  solution  to  the  primal  having  a 
partial  flow  of  I    (<j>)   or  I    (<J>)   in  the  intersected  primal  arc   (y><$). 

The  arcs  whose  lengths  have  been  changed  either  remain  as  members 

of  the  tree  of  shortest  routes  from  node  1  to  all  nodes  (the  early 

tree)  or  become  members  in  the  next  round  of  V.   calculations.  Only 

when  all  are  members  of  an  early  tree  will  the  condition  V.  -  V.  = 

J    1 

V   -  V,   be  satisfied  if  it  was  not  satisfied  initially, 
m    k 

The  forcing  of  the  reverse  arc  of  a  pair  to  a  value  of  -£  (<j>) 
or  -l    (<}>)   does  not  create  a  primal  infeasibility  (that  is,  a  cycle  of 
negative  length  in  the  dual)  for  the  following  reasons.   Suppose  arc 
(i,j)   was  originally  on  the  early  tree.   The  condition  L  .  £  M  . 
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means  that  the  minimum  possible  value  of   £..   will  be  -M  r .   This 

will  occur  only  if  L  »  =  M  .   in  the  primal  problem.   If  arc  (i,j) 

was  originally  in  the  early  tree  then  the  shortest  route  from  node  i 

to  node  j   is  over  arc   (i,j).   Any  other  route  in  the  dual  network 

will  have  a  length  at  least  as  great  as  I . .  =  M  r .   Thus,  even 

ij    YQ 

specifying  I . .    =   -M  .   does  not  create  a  negative  cycle  with  any  of 

the  routes.   The  worst  that  happens  is  that  a  cycle  of  zero  length 

exists  consisting  of  the  arcs   (i,j)   and   (j,i).   If  arc   (i,j)   is 

now  shortened  to  a  value  of   &.(4>)   or   £.,(<(>)   then  it  remains  on  the 

tree  and  I.  .      will  have  a  value  -JL  (<f>)   or   -A,  (d>)   which  is  larger 
ji  1  1 

than  -M  .   so  no  negative  cycle  is  created. 

Next,  suppose  arc   (i,j)   was  not  originally  in  the  early  tree. 
There  must,  therefore,  be  some  other  route  between  nodes   i  and  j 
which  is  shorter  (it  might  be  arc   (j,i)).   Reducing  the  length  of 
(i,j)   to  &..($)   or   £,(<}>)   causes  it  to  become  an  alternate  member 
of  the  early  tree.   Changing  the  length  of  arc   (j,i)   to  a  length  of 
-A,  (<J>)   or  -l   ($)   creates  a  cycle  of  zero  length  with  arc   (i,j). 
If  the  original  shortest  route  between  nodes   i  and  j   was  not  the 
arc   (j,i)   but  rather  some  other  route  directed  from  node   i   to 
node  j   then  its  length  is  at  least   £,  (<j>)   or  I    (<J>)   or  arc   (i,j) 
would  not  be  the  alternative  shortest  route  after  the  new  V.   computa- 
tions.  Therefore,  at  worst  a  cycle  of  zero  length  would  be  created  by 
setting  I    .      at  a  value  of  -I.  (4>)   or  -9.    (<j>) .   If,  instead,  the 
original  shortest  route  was  some  route  directed  from  node  j   to  node 

i  then  a  change  in  I . .      has  no  effect  since  it  is  not  a  member  of  a 

ji 

cycle  containing  the  original  shortest  route.  Reduction  of  I  .  to 
the  value  of  l.  (<|>)  or  H-  (<J>)  will  at  worst  create  a  cycle  of  zero 
length  when  combined  with  this  route. 
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9.   COMMENTS  ABOUT  THE  ALGORITHM 

The  intent  of  the  algorithm  of  section  6  is  to  provide  a  means 
of  solving  dual  networks  for  primal  network  optimal  flows  without  the 
worry  of  whether  the  network  is  source-sink  planar.   All  that  is  required 
is  that  the  primal  network  be  drawn  such  that  no  arc  intersects  any  other 
arc  more  than  once.   While  the  algorithm  is  designed  only  to  solve 
maximal  flow  problems,  the  necessary  modifications  for  solving  the 
minimal  flow  problem  are  minor  and  obvious. 

This  algorithm  may  not  be  as  efficient  an  approach  as  the  maximal 
flow  algorithm  of  Ford  and  Fulkerson  even  if  no  pseudo  nodes  are  needed 
because  the  dual  graph  must  be  constructed.   The  existence  of  pseudo 
nodes  could  reduce  efficiency  further  since  several  iterations  of  the 
algorithm  may  be  required  before  an  optimal  flow  is  found  or  the  primal 
network  is  found  to  be  infeasible.   Careful  construction  of  the  two- 
dimensional  representation  of  the  primal  can,  however,  reduce  pseudo 
nodes  to  a  minimum. 

The  algorithm  is  quite  useful  and  efficient  for  sensitivity 
analyses  and  parametric  studies  of  network  flow  as  a  function  of  primal 
arc  capacities  and  for  problems  involving  capacity  expansion  and  reduc- 
tion (see,  for  example,  Doulliez  and  Rao  (1971)  and  McMasters  and  Mustin 
(1970)).   In  such  problems,  the  dual  network  needs  to  be  constructed 
only  once.  The  changes  in  primal  arc  capacities  then  appear  as  arc 
length  changes  in  the  dual  and  a  new  shortest  route  can  be  quickly 
determined  when  necessary  for  each  change. 
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